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Abstract. Studies of the kinematics of gaseous halos in some spiral galaxies 
have shown that the rotational velocity of the gas decreases with height above 
the midplane (z). This vertical velocity gradient has been observed in both 
the neutral and ionized components of halos. We present imaging Fabry-Perot 
Ha observations of the edge-on galaxy NGC 5775. These observations allow 
us to study the variation in the rotation curve (and, to some extent, the radial 
profile) of the ionized gas component as a function of z. We have used an iterative 
technique to infer the major axis rotation curve from the ionized gas observations 
as well as CO data. A preliminary analysis of position-velocity (PV) diagrams 
parallel to the major axis implies the presence of a vertical velocity gradient. 
The magnitude of this gradie nt is t e ntativ ely f ound to be 
consistent with the results of lRandl l)2000() and lTiillmann et alJ l)200C 



Introduction 



Deep observations of the gaseous halos in some spira l galaxies (e.g., NGC 891, 
ISwaters. Sancish fc van der Hulstl Il997l : UGC 7321. iMatthews fc Woodl fcoOS: 
NGC 2403. iFraternali et alJ l2002 ) have revealed a decrease in rotational veloc- 
ity with height above the midplane (z). This decrease has been observed in both 
th e neutral and diff use ionized gas (DIG) halo components. The fountain model 
of lBregmanl fl980) may explain this behavior because as gas is lifted into the 
halo, it feels a weaker gravitational potential, migrates radially outward, and 
thus its rotation speed drops in order to conserve angular momentum. Thus far, 
the global effects on the kinematics of an individual galaxy have been simulated 
by tr eating the fountain as a purely ballistic pro cess (ICollins. Benjamin. & Rand 
2002]), and as a purely hydrodynamic process ([Beniaminl 2000; Ciotti, this vol- 
ume) . Further detail and comparisons with available data are provided by Rand 
(this volume). 
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While H I observations provide data that can be used to investigate the 
vertical velocity gradient in a spatially complete manner, complementary pro- 
grams in the optical typically make use of slit spectroscopy. Therefore, two- 
dimensional velocity data for the DIG component are usually not available. A 
prime candidate for slit spectroscopy has been NGC 5775, because it is rela- 
tively nearby, nearly edge-on (D = 24.8 Mpc, i = 86°: Ilrwm1ll99^. and shows 
a rem a rkably thick and bright DIG layer with vertical filaments (|Collins et al.l 
2000). iRandl (|2000h found that the velocities measured along two slits parallel 
to the minor ax i s were consistent with a decrease in rotational velocity with z. 
iTullmann et al.1 (I2000T ) found the same trend with a third slit, also parallel to the 
minor axis. In this paper, we present Fabry-Perot observations of the ionized 
gas in NGC 5775. These observations enable us to study the rotation speed as 
a function of z in a manner similar to studies of the vertical velocity gradient 
in the neutral component of galaxy halos. We briefly describe the observations 
and data reductions in §2, describe an analysis of the velocity field in §3, and 
summarize in §4. 

2. Observations and Data Reduction 

Details of the observations and data reduction techniques will be described in 
detail in a forthcoming paper. For the sake of brevity, only the most impor- 
tant details are included here. Data were obtained during the nights of 11-13 
April 2001 at the Anglo- Australian Observatory (AAT), using the TAURUS-II 
Fabry-Perot interferometer (order 379) in conjunction with the MIT/LL 2kx4k 
CCD and an order blocking filter (6601/15). The spectral resolution is FWHM 
~ 0.5 A = 22.9 km s _1 , and the free spectral range (FSR) is 17.4A. Standard 
data reduction steps (bias and flat field corrections and cosmic ray removal) were 
performed using IRAF 1 . Sky-line removal, intensity and wavelength calibrations, 
construction of a data cube, and continuum subtraction were performed using 
scripts written in MATLAB. The final cube was smoothed to a spatial resolu- 
tion of 3.67 arcseconds. The noise in the channel maps was measured to be 
4.93 x 10 -18 erg cm -2 s -1 arcsec -2 (assuming T = 10 4 K, this corresponds to an 
emission measure of 2.47 pc cm -6 ). 

3. Analysis and Modeling 

In order to characterize the kinematics of the DIG component in the halo of NGC 
5775, one can glean some basic tre nds from the velocity field (m oment- 1 map). 
However, as several authors (e.g. iKregel &: van der Kruitll2004h have pointed 
out, this information is of limited use in the study of edge-on galaxies. For a 
given resolution element in the plane of the sky, the observed velocity profile 
contains contributions from every detectable location in the disk intercepted by 
the line of sight (LOS). In an edge-on galaxy, the result is that each velocity 



1 IRAF is distributed by the National Optical Astronomy Observatories, which are operated by 
the Association of Universities for Research in Astronomy, Inc., under cooperative agreement 
with the National Science Foundation 
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profile contains contributions from many different galactocentric radii (R). In 
general, the gas at these radii is characterized by different gas densities and 
different rotational velocities v(R). It is clear, then, that the recovery of a 
rotation curve in an edge-on spiral is more complicated than simply calculating 
the mean velocity. 

Several methods have been suggested for recovering the rotation curve of 
a galaxy. In particular, two useful methods in the study of edge-on spi rals are 
the envelope tracing method ( Sancisi k, Allenlll979l : ISofue k, Rubinl200lh , which 
calculates the rotation curve using the high-v elocity edge of a position- velocity 
(PV) diagram, and the "iteration method" (fTakamiv a k Sofud 12002). which 
automates the procedure of generating a model galaxy (with specified rotation 
curve and density distribution) that best matches the kinematic properties of the 
observed galaxy. The benefit of using the iteration method is a more accurate 
recovery of the rotation curve at small radii, which, because of beam smearing 
and rapidly changing densities and velocities, is typically underestimated by the 
envelope tracing method. 

We have written MATLAB scripts to implement the envelope tracing and 
iteration methods. Space constraints prohibit a description of the algorithms, 
which will be presented in detail in a forthcoming paper. We note, however, that 
optimizing the performance of the iteration method is critically dependent on 
the selection of an appropriate radial density profile. In addition, the signal-to- 
noise and velocity dispersion in the model galaxy must match those of the data. 
These requirements are imposed by the fact that the iteration method seeks to 
generate a model PV diagram which best matches a corresponding PV diagram 
derived from the data, using rotation curves generated from these diagrams by 
the envelope tracing method as a convergence criterion. 

The envelope tracing and iteration methods work under the assumption 
that the disk is transparent. When working in the optical, this assumption is 
not necessarily valid. Indeed, optical images of NGC 5775 reveal a prominent 
dust lane running parallel to the major axis. While the location of the dust lane 
indicates that it may not interfere with calculating the rotation curve along the 
major axis, extinction is still a concern. With that in mind, we have calculated 
a major axis rotation curve using CO 2-1 data ( Lee et ai1l200lh kindly provided 
by S.-W. Lee. We assume that the spatial and velocity distribution of CO and 
DIG is identical in the disk so that we can compare rotation curves derived 
separately from the two individual data sets. 

Figure 1 shows the major axis rotation curve of NGC 5775, derived from 
the CO and Ha data, using the iteration method. In both cases, the points 
were determined by clipping the velocity profiles at 3 times the rms noise in 
the PV diagram. Error bars were determined by clipping the velocity profiles 
at 2 and 4 times the rms noise a nd running the program again. The rotation 
curve determined bv llrwir] (1994) for the neutral gas component is plotted for 
reference. 

The CO data have a much larger beam (21") than the Ha data (3.67"). 
Thus, it should not be too surprising that the rotation curve derived from the 
CO data is much smoother. We only show the Ha rotation curve obtained for the 
approaching side because the receding side rotation curve was found to be very 
low, perhaps due to extinction obscuring the gas on the line of nodes. Further 
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Figure 1. Major axis rotation curve of NGC 5775, derived from (a) CO 
data for the approaching side (crosses) and receding side (circles), a nd (b) 
Ha da ta for the approaching side (triangles). The rotation curve from llrwinl 
(1994) is plotted for reference (dashed line). 



analysis will be required to explain the kinematics of the gas on the receding 
side of the major axis. 

To check for a dropoff in rotational velocity with z, we have modified the 
GIPSY task GALMOD to include a linear vertical velocity gradient: 

v(R,z) = ( «(£* = 0)-S(M-«,) iorz>z ( 
L v(R,z = U) lor z < zq, 

where zq is the scale height of the vertical gas distribution, assumed to be ex- 
ponential, and dv/dz is a constant parameter in the model. We have generated 
galaxy models using this code, and visually compared model PV diagrams from 
the halo with corresponding diagrams obtained from the data. It is important 
to note that in this preliminary study, the models are constructed under the 
assumption that the shape of the radial density profile is constant with z. This 
assumption will not be true if halo gas is, for example, radially outflowing. 

To begin the analysis, we have examined PV diagrams at \z\ = 20 arcsec 
above the midplane, on both the NE and SW sides of the disk. This height was 
selected to eliminate concerns abou t extinction, and to a void confusing the halo 
with the edge of the inclined disk. Collins et al.l (|2000l ) report that extinction 
effects should become negligible at z ~ 600 pc, or 5-6 arcsec. At an inclination 
of 86°, and assuming an optical radius of 120 arcsec, emission from gas at z = 
(in the galaxy frame) should be negligible at projected z-heights greater than 
10 arcsec. 

In Figure 2, we present a sequence of PV diagrams. Note that the receding 
(Offset RA > 0) side will not be well matched, as we have not yet derived a 
reliable rotation curve for that side. In the upper left panel, the major axis PV 
diagram shows the extent to which we were able to match the model (black) 
to the data (white) in the midplane. Next, in the upper right panel, the PV 
diagram from a slit parallel to the major axis, but offset by 20 arcsec to the SW, 
is displayed and compared to the same model. No vertical velocity gradient has 
been added. The inner part of the PV diagram from the data is shallower, and 
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Figure 2. Position- velocity diagrams from the Ha data (white) and galaxy 
models (black). Upper left: Major axis; contours run from 5 to 505xrms, in 
increments of 50xrms. Upper right: z = —20 arcsec; contours run from 5 to 
55xrms, in increments of 5xrms. Lower left: z — —20 arcsec; same contour 
values; v sys has been offset by 20 km s _1 . Lower right: z = —20 arcsec; same 
contour values; v sys has been offset by 10 km s _1 , and a vertical velocity 
gradient of 2 km s _1 arcsec _1 has been added. 

reaches lower velocities than are seen in the model. In the lower left panel, we 
have added 20 km s _1 to the systemic velocity in the model. The inner region is 
much steeper in the model, but this may be explained by a change in the radial 
density profile (the radial profile and rotation curve are very difficult to decouple 
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in the inner region). Therefore, we cannot yet rule out a shift in the systemic 
velocity of the halo. Finally, in the lower right panel, we display a model with 
a vertical velocity gradient of 2 km s -1 arcsec -1 . The systemic velocity has 
been offset by 10 km s -1 to achieve a better match. Again, a modification of 
the shape of the radial density profile may be required. It is not yet clear which 
combination of parameters produces the best model, but this final model appears 
most like the data. Results from the NE side of the major axis are similar in 
appearance. 

4. Summary and Future Work 

The dumpiness of the ionized gas distribution, and the possibility of a change 
in the radial density profile with height, makes it difficult to choose the best 
model. At this stage of our analysis, a vertical velocity gradient of 2 km s -1 
arcsec -1 is consistent with the observed velocity data. However, an offset in the 
systemic velocity of the halo without a changing rotation curve appears to fit 
the data as well. Clearly, further analysis is in order. 

As our analysis continues, we are examining PV diagrams along additional 
slices parallel to the major axis. The extra information from these slices may 
enable us to distinguish between a systemic velocity shift in the halo and the 
changing rotation curve predicted by the fountain model. Moreover, if fountain- 
type gas is flowing radially outward as it is lifted into the halo, the effects may 
be more easily recognized as a change in the radial distribution of gas with 
increasing z. This possibility will be explored. We may also be able to detect a 
lagging halo directly by deriving rotation curves for gas at z > 0. 

Acknowledgments. This material is based on work parially supported by 
the National Science Foundation under Grant No. AST 99-86113. We thank 
Filippo Fraternali for the use of his GIPSY task. We are also grateful to Siow- 
Wang Lee for providing the CO 2-1 data cube, and Judith Irwin for providing 
an H I cube. 

References 

Benjamin, R. A. 2000, Rev. Mexicana Astron. Astrofis. Ser. Conf., 9, 256 
Bregman, J. N. 1980, ApJ, 236, 577 

Collins, J. A., Rand, R. J., Duric, N., & Walterbos, R. A. M. 2000, ApJ, 536, 645 
Collins, J. A., Benjamin, R. A., & Rand, R. J. 2002, ApJ, 578, 98 
Fraternali, F., van Moorsel, G., Sancisi, R., & Oosterloo, T. 2002, AJ, 123, 3124 
Irwin, J. A. 1994, ApJ, 429, 618 

Kregel, M. & van der Kruit, P. C. 2004, MNRAS, 352, 787 

Lee, S.-W., Irwin, J. A., Dettmar, R.-J., Cunningham, C. T., Golla, G., & Wang, Q. D. 

2001, A&A, 377, 759 
Matthews, L. D. & Wood, K. 2003, ApJ, 593, 721 
Rand, R. J. 2000, ApJ, 537, L13 
Sancisi, R. & Allen, R. J. 1979, A&A, 74, 73 
Sofue, Y. & Rubin, V. 2001, ARAA, 39, 137 

Swaters, R. A., Sancisi, R., & van der Hulst, J. M. 1997, ApJ, 491, 140 
Takamiya, T. & Sofuc, Y. 2002, ApJ, 576, L15 

Tiillmann, R., Dettmar, R.-J., Soida, M., Urbanik, M., & Rossa, J. 2000, A&A, 364, 
L36 



